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The  dynamics  of  separation  mitigation  with  asymmetric  dielectric  barrier  discharges  is  explored  by 
considering  the  gas  flow  past  a  flat  plate  at  an  angle  of  attack.  A  self-consistent  model  utilizing 
motion  of  electrons,  ions,  and  neutrals  is  employed  to  couple  the  electric  force  field  to  the 
momentum  of  the  fluid.  The  charge  separation  and  concomitant  electric  field  yield  a  time-averaged 
body  force  which  is  oriented  predominantly  downstream,  with  a  smaller  transverse  component 
towards  the  wall.  This  induces  a  wall-jet-like  feature  that  effectively  eliminates  the  separation 
bubble.  The  impact  of  several  geometric  and  electrical  operating  parameters  is  elucidated.  ©  2006 
American  Institute  of  Physics.  [DOI:  10.1063/1.2187951] 


Active  flow  control  is  vital  to  many  practical  applica¬ 
tions.  Examples  include  low-pressure  turbines  in 
propulsion,1  external  aerodynamics  about  an  aircraft  wing 
where  an  efficient  passive2  or  active  control3  is  important  to 
mitigate  performance  deficiencies  at  off-design  conditions. 
Specifically,  plasma-based  actuators  exhibit  several  potential 
benefits  in  active  flow  control  applications,  including  ab¬ 
sence  of  moving  parts,  rapid  on-off  deployment,  and  attrac¬ 
tive  self-limiting  characteristics.  In  recent  years,  experimen¬ 
tal  observations  have  shown  the  capability  of  dielectric 
barrier  devices,  operating  at  relatively  low  power  levels,  to 
suppress  separation  in  a  wide  range  of  applications  even  at 
atmospheric  pressures.4,5  The  discharge  may  be  characterized 
as  a  transient  microdischarge.  The  time  scale  for  this  dis¬ 
charge  is  smaller  than  microseconds  and  that  for  the  neutral 
gas  is  in  milliseconds.  The  present  authors  have  developed  a 
self-consistent  multifluid  model  for  volume6  and  surface7,8 
barrier  discharges.  The  model  described  the  interaction  be¬ 
tween  the  electric  field  and  the  fluid  mixture  comprised  of 
ions,  electrons,  and  neutrals  under  quiescent  (no  external 
flow)  condition.  In  this  letter,  we  utilize  the  methodology  to 
study  active  separation  control  using  asymmetric  dielectric 
barrier  plasma  actuator.  The  configuration  considered  is  a  flat 
plate  placed  at  an  angle  of  attack  in  an  oncoming  stream  of 
fluid.  The  separated  region  observed  in  the  vicinity  of  the 
leading  edge  of  the  plate  is  then  subjected  to  control  with  a 
suitably  placed  asymmetric  single  dielectric  barrier  plasma 
actuator.  We  solve  the  equations  governing  dynamics  of  elec¬ 
trons,  ions,  and  fluid  to  obtain  spatiotemporal  profiles  of 
electron  density,  ion  density,  electric  potential,  fluid  density, 
and  fluid  velocity. 

A  typical  electrode  arrangement  for  such  actuation  de¬ 
vices  is  shown  in  Fig.  1  where  the  grounded  electrode  is 
surrounded  by  an  insulator  and  a  voltage  fluctuating  at  rf  is 
applied  to  the  electrode  exposed  to  the  gas.  Here,  the  work¬ 
ing  gas  is  helium  at  a  bulk  pressure  and  temperatures  of 
300  Torr  and  300  K,  respectively.  Positive  helium  ions  are 


formed  through  a  process  of  direct  and  stepwise  ionization, 
and  dissociative  recombination.  The  electric  field  E  exerts  a 
net  force  qE  through  the  space  charge  (q)  separated  plasma 
within  the  dielectric  barrier  discharge  (DBD).  The  param¬ 
eters  controlling  such  force  are  the  applied  voltage,  fre¬ 
quency,  dielectric  characteristics,  the  asymmetric  configura¬ 
tion  of  the  electrodes,  and  the  thickness  of  the  exposed 
electrode.  The  effect  of  these  parameters  on  separation  con¬ 
trol  mechanism  is  also  characterized  in  this  letter. 

The  region  simulated  is  10  cm  long  and  5  cm  high.  The 
lower  part  of  the  domain  consists  of  a  0. 1  cm  thick  insulator 
with  dielectric  constant  ed=ae0  while  the  upper  part  is  filled 
with  inert  helium  gas  of  ef=  1.0055  e0,  where  e0  is  permittiv¬ 
ity  of  the  free  space.  The  embedded  electrode  is  grounded 
and  a  sinusoidal  voltage  cf>0  sin(2 irft)  is  applied  to  the 
exposed  electrode.  The  frequency  of  excitation  is  fixed  at 
/=  5  kHz  while  the  influence  of  different  amplitudes,  </>0,  is 
also  explored. 

The  drift-diffusion  form  of  continuity  and  Poisson’s 
equations  for  the  electrons  and  ions  are  solved  with  the  fluid 
momentum  and  continuity  equations  as  described  in  (Ref.  7 
and  8).  The  bulk  density  of  the  helium  is  taken  to  be 
1.79  X  10  4  g/ cm3,  and  the  viscosity  is  assumed  to  be  1.9 
poise.  The  self-consistent  formulation  is  solved  using  a 
variational  formulation  based  on  Galerkin  finite-element 
method9  to  obtain  electron  and  ion  density,  electric  potential, 
neutral  velocity,  and  density. 


^Author  to  whom  correspondence  should  be  addressed;  electronic  mail:  FIG.  1 .  Schematic  of  an  asymmetric  single  dielectric  barrier  plasma  actuator 

sroy@kettering.edu  with  an  incident  gas  flow  angle  a. 
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FIG.  2.  (Color  online)  (a)  The  force 
per  unit  volume  e(w,— wg) E  (dyn/cm3) 
and  currents  components  as  function 
of  cat  at  x=1.75  cm  and  y=0.25  cm. 
(b)  Streamlines  and  vectors  of  gas  ve¬ 
locity  at  the  end  of  1  s,  10  cycles,  30 
cycles,  and  50  cycles. 


For  all  cases  reported  here,  the  rf  electrode  extends  from 
v=0.3  to  1.5  cm  at  y  =  0.1  cm,  the  grounded  electrode  is  from 
v=1.48  to  2.75  cm  at  y  =  0,  with  a  0.02  cm  overlap  between 
electrodes  along  the  v  axis.  The  initial  velocity  components 
u0  (along  the  v  direction)  and  u0  (along  the  y  direction)  are 
1000  and  175  cm/s,  respectively,  corresponding  to  an  angle 
of  attack  ( a )  of  approximately  10°.  The  no-slip  condition  is 
assumed  for  the  gas  neutrals  at  the  dielectric  and  exposed 
electrode  surfaces  and  the  velocity  at  the  left  boundary  of 
upper  domain  is  set  to  the  freestream  condition  at  all  times. 
Homogeneous  Neumann  conditions  are  applied  to  the  upper 
edge  of  the  domain.  For  the  charge  equations,  the  total  cur¬ 
rent  continuity  is  ensured  across  the  dielectric  interface;  i.e., 
at  this  location,  conduction,  convection,  and  displacement 
currents  in  the  gas  are  balanced  with  the  displacement  cur¬ 
rent  in  the  dielectric. 

The  present  study  reports  a  parametric  study  in  which 
applied  rf  voltage  <£0,  initial  plasma  density  n0 ,  dielectric 
constant  a ,  and  rf  electrode  shape  are  varied.  Unless  other¬ 
wise  stated,  ft0=1010/cm3,  a  =  3.5,  cf>0= 400  V,  and  rf  elec¬ 
trode  thickness  is  negligible.  Since  the  gas  is  preionized,  a 
reasonably  small  rf  potential  is  assumed  sufficient  to  main¬ 
tain  the  discharge. 

The  electrons  are  repelled  by  the  insulated  electrode  and 
by  the  exposed  electrode  during  most  of  the  positive  and  the 
negative  part  of  the  rf  voltage,  respectively.  The  local  imbal¬ 
ance  arising  from  the  disparate  mobilities  of  the  ions  and 
electrons  results  in  a  time-varying  charge  separation  distribu¬ 
tion  and  establishment  of  a  consistent  electric  field.  The  gen¬ 
erated  electric  field  restricts  movement  of  electrons  away 
from  the  domain.  At  all  phases  in  the  cycle,  the  peak  mag¬ 
nitude  of  the  force  is  located  downstream  of  the  rf  electrode. 
Figure  2(a)  plots  the  calculated  forces  per  unit  volume 
e{ni~ne) E,  where  the  electric  field  E  =  -V</>,  and  the  dis¬ 
charge  current  as  function  of  cot  at  x=1.75  cm  and 
y  =  0.25  cm.  Although  the  x  component  of  the  force  takes 
both  positive  and  negative  values  during  positive  and  nega¬ 
tive  parts  of  the  cycle,  respectively,  its  magnitude  during  the 
former  phase  is  nearly  three  times  larger  than  in  the  latter. 
Magnitude  of  force  components  is  the  highest  around  posi¬ 
tive  peak  of  the  cycle  and  almost  zero  at  the  end  of  the  cycle. 
The  time  average  of  the  force  at  this  point  is  oriented  in  the 
positive  v  and  negative  y  direction.  Consequently,  collision 

effects  ensure  that  the  dominant  neutral  species  experience 
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an  accelerating  body  force  downstream  and  towards  the  di¬ 
electric  surface.  Since  fluid  cannot  penetrate  the  dielectric,  it 
is  turned  parallel  to  the  surface  of  the  dielectric.  The  net 
effect  as  discussed  further  below  is  therefore  a  gradual  miti¬ 
gation  of  flow  separation.  Ensuring  the  proper  ratio  between 
accelerating  and  retarding  forces  is  the  key  to  successful  ap¬ 
plication  of  DBD  for  actuating  effect  at  higher  neutral  gas 
speeds. 

The  effect  of  DBD  is  highlighted  by  considering  the  flow 
field.  The  initial  field,  without  plasma  actuator,  is  shown  in 
the  first  frame  of  Fig.  2(b),  which  depicts  the  stream  traces  of 
neutral  velocity  at  1  s  (roughly  1000  times  the  characteristic 
flow  time  scale)  after  initiation  of  the  flow.  A  separation 
bubble  develops  over  the  exposed  surface  degrading  perfor¬ 
mance  due  to  increased  drag  and  turbulent  flow  structures. 
When  the  plasma  actuator  is  switched  on,  a  highly  transient 
process  is  initiated  as  shown  in  the  last  three  frames  of  Fig. 
2(b),  which  depicts  stream  traces  of  neutral  velocity  after  10, 
30,  and  50  cycles,  respectively.  The  attachment  process 
downstream  of  the  rf  electrode  progresses  successively:  the 
attachment  point  is  at  6.5  cm  after  10  cycles,  moving  to  8  cm 
after  30  cycles  and  finally  beyond  the  computational  domain 
after  50  cycles.  The  transfer  of  momentum  results  in  a  near¬ 
wall  energized  flow  of  neutrals,  altering  the  dynamics  of  the 
inertial  and  adverse  pressure  gradient  terms  to  eliminate  the 
separation  bubble.  Since  the  neutral  density  is  nearly  107 
times  higher  than  that  of  the  plasma,  therefore,  quick  acting 
controls  may  be  limited  in  high  speed  flows  and  inertia  ef- 


(a) 


(b) 


FIG.  3.  (Color  online)  (a)  Velocity  u  (cm/s)  and  the  force  Fx  (dyn/cm3)  as 
a  function  of  rf  voltage  (f>0  (in  volts),  (b)  Gas  velocity  u  as  a  function  of  x  in 
cm.  Different  lines  are  at  y=0.25  cm  for  initial  plasma  density 
5Xl010/cm3  (Al),  lX10n/cm3  (A2),  dielectric  constant  6=7  (Bl),  and 
6=10.5  (B2). 
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FIG.  4.  (Color  online)  (a)  Streamlines 
and  vectors  of  gas  velocity  for  the 
cases  with  5  X  1010/cm3,  0O= 450  Y, 
finite  thickness  (0.05  cm)  and  arc 
shape  rf  electrode,  (b)  Time-averaged 
force  distribution  showing  the  effect  of 
electrode  radius,  x  and  y  are  in  cm. 


(b) 


fects  are  reflected  in  the  response  time  taken  to  attach  the 
flow  over  full  length  of  the  dielectric. 

The  effect  of  excitation  amplitude  was  explored  by  con¬ 
sidering  </>0  values  of  450,  400,  and  300  V.  For  these  three 
values,  the  correlation  between  the  stream  wise  component  of 
velocity  and  body  force  is  highlighted  in  Fig.  3(a)  which 
plots  them  against  amplitude  of  excitation.  The  magnitudes 
of  the  two  quantities  are  correlated  with  applied  voltage  in 
showing  the  relationship  G~C(f) ™  (cf.  Ref.  5).  For  G  =  u  in 
cm/s  and  </>0in  volts  we  determined  C=77.0  and  m  =  0.4, 
while  for  G  =  FX ,  C=0.03  and  m=  1.2.  This  indicates  that  due 
to  viscous  drag  even  as  the  force  increases,  the  resulting 
imparted  fluid  velocity  does  not  keep  growing  indefinitely 
and  may  reach  a  point  of  diminishing  return  beyond  certain 
optimum  electrical  input.  For  all  cases,  the  positive  peak  of 
the  force  occurs  near  v=  1.75  cm,  i.e.,  just  downstream  of  the 
trailing  edge  of  the  exposed  electrode,  while  negative  peak  is 
observed  at  x~3  cm.  As  increases  by  50%,  the  stream- 
wise  force  increases  by  a  commensurate  amount.  The  mag¬ 
nitudes  of  both  the  components  of  the  force  decrease  sharply 
with  y  and  are  negligible  beyond  x=4  cm. 

The  effect  of  dielectric  constant  and  initial  plasma  den¬ 
sity  is  explored  in  Fig.  3(b),  which  shows  stream  wise  neutral 
velocity  for  initial  plasma  density  5  X  10 10  cm-3  (curve  Al) 
and  1011  cm-3  (A2),  and  dielectric  constant  a  =  7  (Bl)  and 
10.5  (B2),  respectively.  The  higher  value  of  initial  density 
results  in  larger  charge  separation  and  hence  higher  gener¬ 
ated  electrostatic  field.  Consistently,  the  larger  force  e(nt 
-ne) E  yields  more  energetic  plasma  and  neutral  velocities. 
Variation  of  dielectric  constant  does  not  have  a  significant 
effect  on  the  velocities.  A  summary  comparison  of  perfor¬ 
mance  of  a  single  pair  actuator  under  four  different  design 
conditions  is  plotted  in  Fig.  4(a).  Unless  otherwise  stated  for 
these  cases  ?z0  =  5  X  1010/cm3,  <fi0= 450  V.  During  simula¬ 
tions  it  was  found  (not  shown)  that  the  flow  is  smoother  for 
plasma  density  5  X  10 10 /cm3  than  those  for  plasma  densities 
1010/cm3  and  10n/cm3.  The  flow  is  best  attached  to  the 
surface  in  the  case  of  arcshape  rf  electrode  case  which  is  due 
to  increase  in  the  x  component  of  the  force  e{nf*n^ E  as 
observed  earlier.  Figure  4(b)  demonstrates  the  effect  of  elec¬ 
trode  radius  of  curvature.  As  compared  to  the  sharp  corner 
raised  electrode,  the  cycle  averaged  force  magnitude  nearly 
doubles  for  an  arc  shaped  electrode  with  a  small  cusp  at  the 
tip.  This  confirms  the  experimental  observation  reported  in 
Ref.  5.  It  is  important  to  identify  that  in  air,  negative  ions  and 


metastable  species  may  be  present.  However,  the  essential 
characteristics  of  the  electric  body  force  development  and  its 
distribution,  and  the  evolution  of  separation  mitigation  will 
be  quite  similar  to  that  described  here. 

In  conclusion,  the  impact  of  momentum  transfer  and  the 
consequent  enhancement  of  near- wall  neutral  velocity  on  the 
flow  structure  have  been  explored.  A  net  streamwise-oriented 
force  is  generated,  which  acts  on  plasma  predominantly  in 
the  region  downstream  of  the  rf  electrode.  Momentum  trans¬ 
fer  from  charged  particles  to  the  neutrals  results  in  an  ener¬ 
gized  near-wall  flow,  facilitating  the  elimination  of  a  separa¬ 
tion  bubble.  Several  parameters  of  the  actuator  affect 
separation  control  characteristics,  including  the  length  of  the 
rf  electrode,  its  shape  and  location  relative  to  the  embedded 
electrode,  and  properties  of  the  dielectric,  rf  voltage  charac¬ 
teristics  and  initial  plasma  density.  Simulation  confirms  the 
experimental  observation  that  the  force  magnitude  increases 
with  decreasing  radius  of  electrode.  However,  scalability  to 
large  systems  remains  a  concern  since  induced  force  and 
fluid  velocity  do  not  increase  indefinitely. 

In  real  applications,  air  may  be  used  and  water  or  water 
vapor  may  be  present.  This  will  affect  the  chemistry  of  ions 
formation.  The  separation  can  be  controlled  in  either  case 
due  to  the  directional  body  force.  A  thin  insulating  layer  over 
the  upper  electrode  may  be  also  used  to  increase  the  durabil¬ 
ity  of  the  electrode.  Optimization  methodologies  are  neces¬ 
sary  to  generate  such  effective  designs. 

This  work  is  partially  supported  by  Air  Force  Research 
Laboratory  Contract  No.  F33615-98-D-3210  and  support 
from  the  AFOSR  Grant  No.  FA9550-05- 1-0074  (Monitor: 
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